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Abstract. Smoluchowski's coagulation equations can be used as elementary mathemati- 
cal models for the formation of polymers. We review here some recent contributions on a 
variation of this model in which the number of aggregations for each atom is a priori lim- 
ited. Macroscopic results in the deterministic setting can be explained at the microscopic 
level by considering a version of stochastic coalescence with limited aggregations, which 
can be related to the so-called random configuration model of random graph theory. 
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1. Introduction 

In 1916, Smoluchowski introduced an infinite system of differential equations 
which are meant to describe the evolution of the concentrations of particles in a 
medium such that pairs of particles merge as time passes. Since then, this area 
has been intensively studied by physicists, chemists, and mathematicians. See for 
instance the monographs by Drake [S] and by Dubovski [lOj . and the survey by 
Aldous [T]; of course not to mention the numerous references therein and further 
recent works. 

In the original model of Smoluchowski, particles are entirely determined by their 
sizes. Here, we will be interested in the case when they have also arms (or stubs) 
which serve to perform aggregations. In other words, each particle has initially 
a certain number of potential links which are consumed when this particle takes 
part to a coagulation event. The effect is that the total number of aggregations 
involving a given particle is limited by its initial number of arms, and we shall see 
that this has a significant impact on concentrations. 

The purpose of this text is to survey some recent contributions to this topic, 
including explicit resolution and discussion of phenomena such as gelation and self- 
organized criticality. Formulas for the limiting concentrations when time goes to 
infinity hint at connexions with certain branching processes. These are explained 
by the study of the microscopic model. The latter is a close relative to the so-called 
random configuration model, a random multigraph in which the degree sequence 
is pre-described. 

The rest of this paper is organized as follows. Section 2 contains some back- 
ground on Smoluchowski's equations, with a focus on the multiplicative kernel 
and the phenomenon of gelation. Section 3 is devoted to the system of coagula- 
tion equations with limited aggregations; specifically we shall discuss its explicit 
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solutions, another phenomenon of gelation, and the existence of non-degenerate 
limiting concentrations. Section 4 deals with the microscopic model of stochastic 
coalescent with limited aggregations and its relation with random configurations in 
the subcritical case when gelation does not occur. Finally we shall briefly present 
some recent results of Merle and Normand in the super-critical case, in particular 
an interesting phenomenon of self-organized criticality which arises in this setting. 



2. Smoluchowski's coagulation equation in a nutshell 

2.1. General setting. Consider a medium in which particles coalesce pairwise 
as time passes. Typically, each particle has a mass (or size), say to > 0, and when 
two particles with masses m and to' coalesce, they are replaced by a single particle 
with mass to -|- to'. Smoluchoswki proposed a simple mean- field model for the 
evolution of the concentration Ct{m) of particles of mass to at time t >0. Assume 
for simplicity that masses only take integer values. The dynamics are governed by 
a symmetric kernel k : N* x N* — > [0, oo) which specifies the rate at which two 
particles with masses to and to' coalesce. Smoluchowski's coagulation equation 
reads for every to G N* 



m—l oo 

dtCt{m) — — Ct{rn')ct{m — m')K[m\m — m')~Ct{m) Ct(m')K{ni,m'). (1) 

m' — 1 m' — l 

In the right hand side, the first term accounts for the creation of particles of mass 
TO as result of the coagulation of a particle of mass to' < to and a particle of 
mass TO — to' (and the factor 1/2 is due to an obvious symmetry), and the second 
term for the disappearance of particles of mass to after a coagulation with another 
particle. 

It is often convenient to rewrite ([T]) in the equivalent form 

^ oo oo 

5t(Q,/) = -^ Y,{f{m + m')- f{m)- f{m'))K{m,m')ct{m)ct{m'), (2) 

m— 1 m' — l 

where / : N* — M has finite support and (ct, /) — X]m=i /(™)ct('Ti). The informal 
description of the model might suggest that the average mass of particles in the 
medium 

oo 

(ct. Id) = ^ mct{m) 

m— 1 

should be conserved as time passes, and indeed if we could take / = Id in ([2]), then 
we would get 9t(ct,Id) = 0. However this is not always true; more precisely the 
gelation time 

Tgei-inf{t>0:(ct,ld)^(co,ld)} 

is known to be finite for a large class of coagulation kernels (see e.g. [131 E]). 

Informally, the phenomenon of gelation should be interpreted as the formation 
of giant particles which are not taken into account in the average mass of particles 
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in the medium. The prototype of kernels for which this occurs is the multiplicative 
one; it is discussed next. 

2.2. The multiplicative kernel and gelation phenomenon. We suppose 
here that K.{m,m') = m ■ m', and further, with little loss of generality, that the 
average mass of particles in the media at the initial time i = is (cq. Id) = 1. Then 
for t < Tgci, Smoluchowski equation ([1} takes the simpler form 



^ m — 1 

dtCt{m) — — Ct{m')ct{m — m')m'{m — m') — mct{m). (3) 

m' — 1 

It is easy to solve this system step by step. Indeed, consider for simplicity mono- 
disperse initial conditions, viz. co(m) — lm=i- Then for t < Tgd and m — 1, 
one finds dct(l) — — ct(l)di and thus Cf(l) = e~*. Next for m = 2 one gets by 
substitution 

dQ(2)= Qe-2*_2ct(2)^ dt 
and then Ct(2) = ^ie^^*. By induction, one arrives at the solution due to McLeod 

Q(m) ^t^-^m^-^e-"*/™!, meN*. (4) 

We stress that this is only valid before the gelation time Tgd, which is a finite 
quantity. Indeed {ct , Id) is then given by 

to! 

m=l 

and one can check that 9{t) = 1 for i < 1 whilcQ < 9{t) < 1 for i > 1. In 
particular, Tgd — 1. 

Considering generating functions provides an alternative way for solving Equa- 
tions for f < Tgoi; see e.g. [^8 and the references therein. Indeed, we can 
interpret as a convolution equation 

^dtpt{m) = *Pt{m) -pt{m) 

where pt{m) — mct{m) is the mass- function of some probability distribution on 
N*. In terms of the generating function gt{x) = J2m=i ^^"'Pti'm), x e [0, 1], this 
yields 

dt9t{x) = -^da^gtix) - xdj,gt{x) = x{gt{x) - l)dxgt{x) . 



^More precisely, observe that 

te{t) = — cxp(-(t - lnt)m) , 

m — l 

SO that ift'<l<t are such that t - Int = t' - \nt' , then te{t) =t'e{t') = t' . This entails that 
e(t) = t'/t is the unique solution in (0, 1) to the equation e't^t*'"^' = e{t). 
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This is a quasi-linear PDE which can be solved by the method of characteristics. 
One finds that for < i < 1 (assuming again mono-disperse initial conditions for 
simplicity), gt : [0, 1] ^> [0, 1] is the inverse bijection oi x i-^ xexp{t{l — x)). It is 
known that this is the generating function of the Borel distribution with parameter 
t, so inverting generating functions gives 

Pt{m) = ^-^ , togN*, 

m! 

and one recovers the solution (g]) of McLeod. 

We conclude this section by discussing briefly solutions after the gelation time. 
In the mono-disperse case co(m) — Um^i, Kokholm [19] established global exis- 
tence and uniqueness of the solution, which is given for t > 1 hy 

cAra) = 1 , TO G N*. 

m\t 

Observe that this does not solve the simpler equation ([3]) and that for t > 1, 
(ct,Id) = 1/t > 9{t). For poly-disperse initial conditions, the gelation time can be 
determined explicitly as 

oo 

Tgci = l/(co, m^) = 1/ X! ™^co(m) ; (6) 

?n — 1 

note that gelation may occur immediately, namely Tgd = whenever the second 
moment of cq is infinite. Existence and uniqueness of the solution on the whole time 
interval [0, oo) have been established rigorously only quite recently by Normand 
and Zambotti [27], although of course several partial results appeared previously 
in the literature; see the comments and references in |27j . 

2.3. Connexion with the Erdos-Renyi random graph model. We then 
study coagulation through a microscope, in the sense that we now consider a mi- 
croscopic stochastic model which follows related dynamics at the infinitesimal scale, 
and whose macroscopic evolution is close to Smoluchowski's coagulation equation. 
More precisely, Marcus 21 and Lushnikov have introduced stochastic models 
for the evolution of finite systems of particles with total mass n, in which any pair 
of particles, say (to, to') coagulate at a rate k(to, to')/?!, independently of the other 
pairs. These are known as Marcus- Lushnikov coalescents; see the survey by Aldous 
[T] for background and further references. Roughly speaking, the connexion be- 
tween the deterministic model of Smoluchowski and the stochastic one of Marcus 
and Lushnikov is that solutions to Smoluchowski's coagulation equations can be 
obtained as hydrodynamic limits of Marcus-Lushnikov processes, see e.g. Cepeda 
and Fournier [7], Eibeck and Wagner |T2], Jeon [TH], Norris [MK^n], etc. 

In the case of the multiplicative kernel k(to, to') = to-to', the Marcus-Lushnikov 
coalescent is closely related to a notorious model of random graph theory which 
was studied first by Erdos and Renyi; see e.g. [Q or [T7]. Fix a number n 3> 1 
of vertices and consider the random graph process such that an edge connecting a 
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given pair of vertices appears at an exponentially distributed time with parameter 
1/n (i.e. with expectation n), independently of the other pairs of vertices. In 
particular, the rate at which an edge appears between a connected component of 
size m and another connected component of size m' is (m • m')/n (this is 1/n times 
the number of edges which can possibly connect these two clusters). In other words, 
the process which records the sizes of the clusters in the random graph process is 
Markovian; more precisely it is a version of the Marcus-Lushnikov multiplicative 
coalescent started from a mono-disperse initial condition. 

Erdos and Renyi pointed at a nowadays classical phase-transition for the size 
of the connected component in the random graph which contains a typical vertex. 
Roughly speaking, the probability that this size remains bounded when n — oo 
is ~ 1 for t < 1, whereas for t > \, this probability is only ^ 9{t) < 1, where 
9(t) is the unique solution in (0, 1) to the equation e*^^~^^ — x. This provides 
a probabilistic explanation of the gelation phenomenon; recall in particular the 
discussion near ([5]). 

More precisely, the random graph model also enables us to recover the solution 
of McLeod before gelation. Informally, it is easily shown that for n ^ 1, the cluster 
of the random graph process at time t < 1 which contains a given vertex resembles 
a Galton- Watson tree with reproduction law given by the Poisson distribution 
with parameter t. In particular the law of its size converges when n — >■ cx) to the 
Borel distribution with parameter t. Roughly speaking, this entails that m times 
the concentration of clusters of size m (because a given cluster of size m appears 
exactly m times as the cluster containing a certain vertex), is close to the mass 
function of the Borel(i) distribution evaluated at m, i.e. i'"~^m"'~^e~™*/™'- This 
is the formula (U) obtained by McLeod. 

We stress that this connection between the random graph process and Smolu- 
chowski's coagulation equations for the multiplicative kernel ceases after time 1, 
i.e. after the phase transition or the gelation time. The reason is intuitively clear: 
in the random graph process, the giant cluster which appears at time 1 interacts 
thereafter with microscopic clusters (informally, more and more microscopic clus- 
ters are absorbed by the giant component), whereas in Smoluchowski's model, the 
gel (i.e. giant components) does not interact with particles. As a matter of fact, 
the random graph process is even more closely related to a variation of Smolu- 
chowski's coagulation equation, namely the Flory's equation, in which interaction 
with the gel is taken into account; see Fournier and Laurengot |16) . 



3. Coagulation equations with limited aggregations 

3.1. Setting. We now describe the modified model with limited aggregations; it 
may be useful for the intuition to think of the formation of polymers as the result 
of the creation of bonds connecting different atoms. The name particle now refers 
to a polymer which consist of atoms (possibly a single one) connected by bonds. 
Each atom has a certain number of potential connexions, which we call arms, and 
which can be paired with arms of other atoms to create bonds. In this setting, the 
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size m of a particle is given by the number of atoms that constitute this polymer, 
and its number of arms is the total number of potential connections which are still 
available for those atoms. 

We write Ct(a, m) for the concentration at time t of particles with size m and a 
arms in the medium, and assume for simplicity that at the initial time, there are 
no bonds connecting atoms, i.e. Ct{a,m) = for all m > 2. We then define 

H{a) = co(a, 1) , a e N* , 

and assume that the first two moments are finite, i.e. 

oo oo 

Aj -.^^a^coia,!) = ^a^ H{a) < oo for j = 1,2. (7) 

a— 1 a—1 

We imagine that when two particles merge, one arm in each particle is used to 
perform a bond, an event which can be schematically represented as 

(a, m) & (a', m') k> (a + a' — 2, m + m) . 

We further assume that the rate at which such aggregation occurs is a ■ a', i.e. 
it coincides with the numbers of possible connexions between these two particles. 
This means that the equations governing the evolution of concentrations are now 

^ a+l m—1 

dtCt{a, m) = — a' Ct{0'' , Tn') • (a — a' + 2)ct{a — a' + 2, to — to') 

a' — 1 m' — l 

oo oo 

- act{a,m) ■ a'ct{a',m'), (8) 

a' — 1 m' — l 

where here again the first term in the right-hand side accounts for the creation of a 
particle (a, to) as the result of the aggregation of a particle (a', to') with a particle 
(a — a' + 2, TO — to'), and the second for the disappearance of a particle (a, to) after 
an aggregation event. 

Particles with no arms have a special role. They only arise from the aggregation 
of pairs of particles each with a single arm, and cannot be involved into further 
coagulation events. The concentration Ct(0, m) of such particles thus increases with 
time and possesses a limit as t — ?> oo. 

3.2. Solution and gelation. Equation ([SJ somewhat resembles ([1]) for the mul- 
tiplicative kernel, and even though it looks more involved, it can also be solved 
explicitly to a large extent, see [5] and 27 for details. 

Recall from ([6|) that the gelation time in Smoluchowski's coagulation equations 
with the multiplicative kernel can be expressed as the inverse of the second moment 
of the initial condition; in particular gelation always occurs. For coagulation with 
limited aggregations, the gelation time is given by 

r oo if ^2 < 2Ai, 
\i/(A2-2Ai) ifA2>2Ai, 
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where Ai and A2 are defined in ([7|). We stress that T can be infinite, a situation 
which does not occur when aggregations are unlimited. The main result of [3] 
solves (|S]) before gelation; we refer to [57] for more complete results including 
post-gelation solutions. The probability measure ly on N defined by 

^1 

plays a crucial role. 

Theorem 3.1. The system ([5]) has a unique solution on [0,T) started from 

co(a,m) = l{™=i}^(a) , a,meN* 
which is given for a, m > 1 by 

ct{a,m) = (« + - + yli t) 1)^^*"' (a + m - 2) , 

aim! 

where v™ denotes the m-th convolution product of the probability measure v. Fur- 
ther 

Q(0,m) = -(l + l/(Aii))i-'"jy*'"(m-2) form>2. 

m(m — 1) 



Let us now briefiy sketch the main lines of the proof of Theorem 13. II Without 
loss of generality, we may suppose that concentrations have been normalized so 
that Ai = I. One starts by showing that for t < T, the mean number of arms at 
time t is simply given by 



00 00 

a—1 m— 1 



and then ([S]) can be re-written as 

a+l m— 1 



OtCt(a, m) = - a Ct(a ,m )[a — a + 2)ct[a — a +2,to — m) Y~\~t — " 



a' — l m' — l 



Next one uses similar techniques of generating functions as in Section 2. That is 
one defines for x, y G [0, 1] 



9t{x,y) = J2Y1 x°'y'''ct{a,m) 



a—l m—1 

and observes that ([5]) then implies that the function kt{x,y) :— dxgt{x,y) solves 
for t < T the quasi-linear PDE 

dtkt{x,y) = (kt{x,y) - j dxkt{x,y) - Y^kt{x,y) . 
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The latter can be solved by the method of characteristics (see e.g. [TS] for back- 
ground), and one gets 

kt{x,y) = {l+t)-^ko{£t{x,y),y)^t-ht{x,y) 



t^ + t' 



where (•, y) : [0, 1] — > [0, 1] is the inverse to the function a; M> (1 + t)x — tko{x, y). 
All what is needed now is to invert the generating function, which is realized thanks 
to Lagrange inversion formula (see, for instance. Section 5.1 in [32J ) , and we arrive 
at the solution of the statement. 

It is interesting to point out that similar arguments also work in a more so- 
phisticated model in which arms have genders (male or female) and only couple of 
arms with different genders can be paired; see '25 . 

3.3. Limiting concentrations. One remarkable feature of the model with lim- 
ited aggregations is that concentrations have a non-trivial limit as time tends to 
oo. In the case without gelation (i.e. T = oo, or equivalently A2 < 2Ai), this 
is immediately seen from Theorem 13.11 The case where gelation occurs is signifi- 
cantly harder, but somewhat surprisingly, one obtains a very similar formula. The 
following statement rephrases Corollary 6.2 in Normand and Zambotti [27] (and 
also corrects a little misprint there); recall that v*™ has been defined in Theorem 

EH 



Theorem 3.2. The concentration Ct{a, m) of particles (a, m) at time t with a G N 
and m > 2, has a limit Coo(a, fn) as t 00. 

(i) When there is no gelation (i.e. T — 00), the latter is given by 

Coo(a,m) = l{a=o} .. v*"\m - 2) . 

m[m — 1) 

(ii) When gelation occurs (i.e. T < 00), let rj denote the unique solution to the 
equation r]gl{7]) = guiv) where gv{x) = '^x'^v{i) is the generating function of v. 
Then one has 

Coo(a, m) = l{a^o} ^. /?"-V*'"(m - 2) , 

^ ' m(m — 1) 

with (3 = r]/g^{Tj) = 1/gUv) > 1- 

It is interesting to mention further that when gelation does not occur, we have 



mcoo(0,m) = ^co(a, 1) 



which means that the total mass of polymers in the terminal state coincides with 
the initial total mass of atoms (which is stronger than absence of gelation) . 

The explicit formula for the limiting concentration has a strong probabilistic 
flavor. More precisely, consider a Galton- Watson process with reproduction law v, 
i.e. a population model in which individuals reproduce independently one from the 
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other, and with the same offspring distribution v. A classical result due to Dwass 
[TT] claims that the law of the size of the total population in such a Galton- Watson 
process with two ancestors is given by 

— i/*"(m-2), m>2. 

TO 

The purpose of the next section is to provide a probabilistic explanation of the 
appearance of Galton- Watson processes with two ancestors in this setting. 



4. Stochastic coalescence with hmited aggregations 

4.1. The random configuration model and giant components. The deep 
link between Smoluchowski's coagulation equations and Marcus-Lushnikov stochas- 
tic coalescent invites us to consider the following stochastic evolution for finite 
systems of polymers. Recall that a particle (polymer) is characterized by its size 
TO and its number a of arms; a pair of particles, say {a,m) and (a',m'), merge 
at rate a • a' to produce a single particle (a + a' — 2,m + m'), independently of 
the other pairs. Assume henceforth that the initial state of the system is purely 
atomic, i.e. all the particles have unit size, and recall that we are mainly inter- 
ested in the terminal state of the system (i.e. in the limiting concentrations for 
the macroscopic model) . We stress that arms of a polymer are only used to create 
bonds with other polymers; two arms in the same polymer do not joint to form a 
bond. In other words, such polymers have a tree structure. 

The evolution of the system of polymers described above is closely related 
to the random configuration model. The latter is a simple stochastic algorithm 
introduced independently by Bollobas 5 and Wormald 33! in the late 70's, which 
aims at producing a random graph on a set of vertices with pre-described degrees; 
see Chapter 7 in [T7] for a detailed account. Typically, we start with some finite 
set V of vertices and assign to each vertex v G V a, degree d{v) G N*. We imagine 
that there are d{v) stubs attached to the vertex v, and join uniformly at random 
pairs of stubs to create edges between those vertices. The resulting multi-graph 
is known as the random configuration model] in general it is not simple, in the 
sense that there may exist loops and multiple edges. However, it is known that, 
loosely speaking, if the degrees d{v) are sufficiently small, then only few clusters 
in a random configuration have loops, multiple edges, or even cycles. If we discard 
this pathology, the random algorithm which produces the random configuration is 
then remarkably close to the stochastic coalescence with limited aggregations. 

The purpose of this section is to show that, at least in appropriate regimes for 
which, informally, the number of arms is not too large, results on random config- 
urations provide a probabilistic explanation of properties of the solutions to the 
deterministic coagulation equations with limited aggregations. In this direction, 
we first recall that MoUoy and Reed [23j i24j have pointed at a phase-transition for 
large sparse random configurations which generalizes that of Erdos and Renyi for 
large random graphs. Specifically, consider for every integer rt > 1 a set Vn of n 
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vertices and a function c?„ : ^ ^ N* that specifies the number of stubs appended 
to each vertex. Introduce the empirical distribution of the number of stubs 

-.^ -#{v e Vn : dn{v) ^ i} , ieN*. 
n 

Assume that the hmit 

hm fin{i) := fi{i) (9) 

n— >oo 

exists for every z > 1, and further that 

oo oo 

hm y i^n{i) = / < oo . (10) 

71— >-00 ^ ^ ^ ^ 

Informahy, Mohoy and Reed have shown that if 

oo 

5]*(*-2VW<0, (11) 

i=l 

then with high probabihty when n — > oo, the size of a typical connected component 
in a random configuration on Vn remains bounded, whereas when (jlip fails, then 
with high probability there exists a giant cluster of size of order n. Observe that 
([TT|) can be rewritten as A2 < 2Ai in the notation ([7]); so we recover the necessary 
and sufficient condition for the absence of gelation. Deeper connexions will be 
pointed at in the next section. 

4.2. The structure of typical polymers in the sub-critical case. In this 
section, we assume that (|9]), (jlOp and ([TT]) holds, so typical clusters in the random 
configuration model remain microscopic as the number of vertices tends to infinity, 
and gelation does not occur. Our purpose is to show that one can recover the 
formula of Theorem 13.21 for the limiting concentration, by statistical analysis of 
those microscopic clusters. This was achieved in here we shall merely provide 
a rather informal sketch of the argument, refereeing to that article for details. 

Consider the random configuration model on a set Vn of n vertices. For every 
oriented edge e of this random graph, we consider the connected component Ce that 
contains this edge, viewed as a combinatorial structure on a set of unlabeled vertices 
(in the sense that two isomorphic subgraphs are identified) with a distinguished 
oriented edge. We are interested in the empirical measure 

en(r) = ^#{e : = r} 

where Sn is the total number of arms (or equivalently Sn is the number of oriented 
edges) and F is a generic combinatorial structure with a distinguished oriented 
edge. 

Recall that v denotes the probability measure given by 



i^ii) = ii + + 1)/Ai , 



i 6 N. 
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The condition (|lll) is equivalent to requiring that v is critical or sub-critical in 
the sense of branching processes, that is ^ 1- Now consider a Galton- 

Watson process started from two ancestors, with reproduction law v. We know 
that the process becomes extinct a.s., and we can consider the two genealogical 
trees of these two ancestors, which we further connect by an oriented edge from 
the first ancestor to the second one. This gives a combinatorial structure with one 
distinguished oriented edge; we write GWi,^2 for its law. 
The main result of [3] is the following: 

Theorem 4.1. Suppose (ITlT) holds. Then 



where the space of probability measures on the discrete set of combinatorial struc- 
tures with one distinguished edge is endowed with the total variation distance. 

Let us briefly present the idea of the proof, which already appears in different 
forms in the literature on random configurations. From the construction of the 
configuration model, we see that sampling an oriented edge uniformly at random 
amounts to picking two arms uniformly at random at the initial stage of the con- 
struction (i.e. before any bond between vertices has been formed) and joining 
them. Let v and v' be the vertices to which these two arms are attached; possibly 
V — v' , but this event is quite unlikely when n ^ 1 and will thus be discarded. 
Note the degree of v, dn{v) does not have the law /x„ normalized to be a probability 
measure, but rather the size-biased version of the latter. It follows from our as- 
sumptions that dn{v) — 1, which represents the number of arms of v which remain 
available after v and v' get connected, converges in distribution to i^, and by an 
obvious symmetry, dn{v') — 1 also converges in distribution to v. More precisely, 
dn{v) — 1 and dn{v') — 1 are asymptotically independent, in the sense that the pair 
{dn{v) — l,dn{v') — 1) converges in law to a pair of independent variables with the 
same law v. 

Then let us work conditionally on dn{v) — 1 = fc and dn{v') ~ 1 — £. In- 
formally, each available arm of v and of v' will then be connected to different 
vertices, say vi,...,Vk and v[,...,v'^ such that, when n oo, the variables 
dnivi) — 1, • • • , dn{vk) — 1, dn{vi) — 1, • • • , dn(v[) — 1 becomc asymptotically inde- 
pendent and have all the same distribution v. By iteration, we see that if we pick 
an oriented edge e uniformly at random, the cluster Ce converges in distribution 
to a random combinatorial structure with law GWiy.2. 

The proof is then completed with an argument of propagation of chaos. Namely, 
we have to check that if we now pick two oriented edges uniformly at random, say 
e and e', then the combinatorial structures Ce and Ce' become asymptotically 
independent when n — > oo. Indeed, it is well-known that propagation of chaos is 
equivalent to the weak convergence of the empirical distributions; see [?Tj . 

We have now all the ingredients needed for a probabilistic (or microscopic) 
explanation of Theorem 13. 2f i). First, the total number of arms in the system is 



lim e„ = GW,,^2 



in probability, 



n 



oo 



oo 



oo 




1 



1 
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Then, roughly speaking, we know from Theorem 14.11 that most clusters in the 
random configuration model are trees; in particular a tree-cluster C of size |C| = ni 
(number of vertices) has exactly m — 1 edges, which in turn correspond to 2 (to — 1) 
arms. If we think of the random configuration as the terminal state of the polymer 
system, then the concentration of polymers of size to is approximatively 

coo(m)^ ^" e„({C:|q^TO}) 
2(to — l)n 

Recall also that thanks to Dwass' formula, 

GW,,2(|C| = to) = -i/*™(TO - 2) , 
to 



so combining with Theorem 14. 1[ we arrive at 



Coo(to) ~ — -— !— — j/*'"(to - 2) , 
to(to — 1) 



which is the formula found in Theorem 



4.3. Formation of gel and self-organized criticality. We now conclude this 
text by presenting succinctly and informally some resent results obtained by Merle 
and Normand (see Chapter 4 in 26 ) in the supercritical case. 

Roughly speaking, stochastic coalescence can be related to Smoluchowski's co- 
agulation equations only up to the time when a giant component appears, which 
corresponds to the gelation time. This stems from the fact that in the stochastic 
framework, giant components interact with smaller components, whereas this is 
not the case in Smoluchowski's setting. Fournier and Laurengot |16j pointed at 
a modified version of stochastic coalescence which circumvents this problem. De- 
noting by M the total mass of particles in the stochastic system, they introduce a 
threshold a{M) with 1 ^ a(M) <^ M, such that as soon as a particle with size 
greater than a{M) appears, it falls into the gel and ceases to interact with smaller 
particles. In other words, particles with size greater than a{M) are seen as giant, 
and we may imagine that they are then removed from the system. Thus only par- 
ticles with mass less than a{M) are allowed to coagulate, and an hydrodynamic 
limit procedure now enables one to relate this modified stochastic coalescence to 
Smoluchowski's coagulation equations for all times. Note also this modified model 
coincides with the original one of Marcus and Lushnikov before the appearance of 
giant particles. 

Merle and Normand exploited this idea for stochastic coalescence with limited 
aggregations. Strictly speaking, they consider a slightly different model in which 
time runs at a different speed, however this is only a minor change made for 
convenience; in particular this does not alter the terminal state of the system. 
So imagine that we start with n atoms with arms. As time passes, polymers are 
formed by the creation of bonds connecting atoms just as in Section 4.1, except 
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that now polymers with size greater than the thresholcJl a{n) fall instantaneously 
into the gel. 

Merle and Normand focus on the total concentration, say m", of atoms (not 
polymers) in solution at time t, i.e. which have not fallen into the gel before time 
t, and on the empirical measure tt" of used arms in the solution. That is nm" 
represents the total number of atoms in solution at time t, and the number of 
atoms in solution at time t for which exactly k arms have already been used is 
nm" X 7r"(fc). They point at the following simple description of the distribution 
of the configuration in solution: conditionally on n" and tt", the configuration in 
solution has the same law as a random configuration model started with n" atoms 
and empirical distribution of arms tt" , and conditioned on having no cluster of size 
greater than a{n). 

They prove that under hypotheses a bit more restrictive than © and PU)) . to" 
converges as n ^ oo to some deterministic nit G [0, 1] and tt" converges to some 
deterministic probability measure on N*. Both rrit and ttj can be characterized 
explicitly in terms of the limiting distribution of arms /i. It should not come as a 
suprise that before the gelation time, Wt is subcritical in the sense of the random 
configuration model, i.e. 

oo 

^z(i-2)^t(z) < 

i=l 

(recall the condition ([TT|) of Molloy and Reed) ; but a much more interesting feature 
is that TTt remains exactly critical after the gelation, i.e. 

oo 

Y^l{i-2)TTt{i)=0. 

i=l 

This provides a nice illustration of self-organized criticality, a general phenomenon 
which is often observed in Statistical Physics (see |2]), but has been rarely be 
proven mathematically. 

Finally, Merle and Normand consider the limits TOqo — hnif_>.oo fnt and tTqc = 
limt_>oo TTj. Specifically, they obtain that 

oo 

TOoo = ^lffJ,{i) , 
i=l 

where r] has been defined in Theorem 13. 2f ii). and that 

moo 

This enables them to repeat the argument at the end of Section 4.2 and to provide 
a probabilistic explanation for the simple formula of Theorem I3.2f ii) which gives 
the limiting concentrations in the supercritical case. 

■^Merle and Normand assume that n^/'^ ^ a(n) <g n, which is slightly less general than the 
condition of Fournier and Laurengot. 
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